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SUMMARY 

A detailed analysis has been made to take account of the 
high cooling -air velocity occurring in high -altitude radiators. 
Methods are developed for determining the heat -transfer rate, 
the pressure drop, and the drag power . Some effects of Mach 
number are shown. Radiator performance charts based on the anal- 
ysis are presented for a vide range of the design variables. The 
application of the charts is shown by an example. 

The performance charts show that the heat -transfer rate for 
a given total -pressure loss is not greatly affected by the high 
airplane velocities but that the necessary total -pressure loss 
and the resulting drag are both greatly increased at high alti- 
tudes. 



INTRODUCTION 

Extensive literature is available relative to the perform- 
ance of ethylene -glycol radiators in the normal range of oper- 
ating conditions. At high altitudes, however, certain effects 
that normally receive but little consideration acquire increased 
importance as a result of the high velocities of the cooling air 
through the tubes. The purpose of this paper is to describe and 
evaluate these effects. 

The theories are outlined on which are based the calcula- 
tions of heat transfer, of friction pressure drop> and of accel- 
eration pressure drop at high Mach numbers. A general differential 
equation for the pressure drop and some approximate solutions 
for the equation are given. Radiator design charts based on the 
simplest of these approximate solutions are included. These charts 
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show, for a range of altitude up to 50,000 feet, the variation 
of heat dissipation per unit frontal area with pressure drop as 
well as the corresponding drag power. An example showing the 
use of the charts is also included. 

This paper is the third of the series on high-altitude cool- 
ing, (See reference 1.) 



SYMBOLS 

V airspeed, feet per second 

V mean flow velocity in radiator tube, feet per second 
L radiator-tube length, feet 

D radiator-tube diameter, feet 

A cross-sectional area of radiator tube, square feet 

x distance along radiator tube, feet 

y distance from axis of radiator tube, feet 

T absolute temperature, °F + 460 

T v absolute temperature of inner surface of radiator tube, 

°F 4- 460 

T q absolute stagnation temperature, 

v 2 



F + 460, T + 



2Jc 



0 * Z r - T * 
w s 



p density of cooling air, elug3 per cubic foot 

B gas constant for air (53.3 X 32.2 Btu per slug per °F) 

p absolute static pressure, pounds per square foot 

q dynamic pressure, pounds per square foot 
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§ c impact pressure, pounds per square fbot 
H total pressure, pounds per square foot 

Ap^ static-pressure change due to shearing force at tube wall, 
pounds per square foot 

Ap m static -pressure change due to momentum change of cooling 
air, pounds per square foot 

y isentropic-expansion exponent (1.4) 

k thermal conductivity, Btu per square foot per second per 
(°F/foot) 

c specific heat of air at constant pressure (7.75 Btu par slug 
P per °F) 

J mechanical equivalent of heat (778 foot -pounds per Btu) 

h heat -transfer coefficient, Btu per second per square foot 
per F 

Q heat -transfer rate, Btu per second 

P Q heat -transfer rate, horsepower per unit frontal area of 
radiator 

drag power, horsepower per unit frontal area of radiator 

a relative density of atmosphere (p/o. 002378) 

|i absolute viscosity of cooling air, slugs per foot per 
second 

R Reynolds number of flow in radiator tube (V p D/u) 

a speed of sound at temperature T, feet per second 

Subscripts: 

o in free stream 

2 just inside radiator-tube entrance 
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just inside radiator -tube exit 



4 



just outside radiator-tube exit after sudden expansion 
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in cooling air after air returns to free -stream static 
pressure 



a 



stagnation 



ANALYSIS OF THE COOLING-AIR FLOW 



The path of the cooling air through the airplane is con- 
veniently considered in four parts. Between the- free stream and 
the inlet to the radiator tubes the air undergoes a change, usu- 
ally an increase, in static pressure and some loss of total pres- 
sure. In flowing through the hot tube the air undergoes a tem- 
perature rise and a total -pressure drop. In passing out cf the 
radiator tubes the air suffers a further loss of total pressure. 
Finally, as the air passes out of the duct into the free stream, 
the static pressure returns to free-stream static pressure at sub- 
stantially constant total pressure. A detailed analysis of the 
processes occurring in these steps follows. 



The temperature at the radiator -tube entrance is given by 
a statement of the law of the conservation of energy for the 
adiabatic flow of a perfect gas 



where = 0.832 x UCT slug per, F per foot-pound 

?.Jc p 

The pressure just inside the radiator-tube entrance, unlike 
the temperature, is not uniquely determined by the airplane speed 
and the tube -entrance velocity. If the flow from the free stream 



Adiabatic Compression 
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up to this point were isentropic, the pressure pg would he given 

by 




The full isentropic pressure rise is generally not developed 
because of skin friction and flaw separation. An estimate of 
the pressure loss must be made for purposes of calculation. 



Heat Transfer 

The usual equation for the heat transfer from a radiator 
tube is derived by equating the heat lost from an elementary 
length of tube wall to the heat gained by the fluid in flowing 
through the element. The heat lost is 

dQ = h(T - T)*D dx (3) 

* W s 

where T is the temperature at any distance x along the 
radiator tube. The expression for the heat-transfer coeffi- 
cient from reference 2, with the value of the constant from 
reference 3, is 



h . 0.0247c p ^ ' (pV) 0 ' 8 (4) 



The heat gained is given by 



2 ' • • 

dQ = —frVCp dT (5) 



If any change in the kinetic energy within the element is neg- 
lected, an equation relating T and x may be obtained by 
eliminating h and dQ among equations (3), (4), and (5), 
If the variation in the fifth root of the viscosity is neglected, 
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the equation nay be integrated ior the initial conditions T = T 2 
when x = 0; thus 

v -0.2 X 
-4(0,0247)R | 

T - T = (T„ - To)e ( 6 ) 
The quantity of heat transferred per unit time, up to point x is 

Q - c p (T , T 2 ) (7) 

Combining equations (6) and (7) results in the usual equation 
for the heat transfer of radiators 



1 - e 



-4(0.0247)R"'°' 2 ^ 



(8) 



The empirical formula for the heat -transfer coefficient given 
hy equation (4) is based on tests at low airspeeds for which 
the heat developed in the laminar sublayer by viscous shearing forces 
is small when compared with the total heat transfer. Because the 
heat generated by viscous shearing forces is quadratically depend- 
ent on the cooling-air velocity, this heat become appreciable at 
high airspeeds. For such cases Crocco given the following rule 

in re f ace 4: "The transference of heat between any object and 

a fluid 3 lowing by it, when it is no longer possible to disre- 
gard the heat developed by friction, is governed by the same law 
that applied when it is negligible, provided the temperature to 
which the fluid is brought by adiabatic arrest is considered as 
its temperature." In equation (3), that is, the actual tempera- 
ture should be replaced by the stagnation temperature 

dQ » h(T w - T S )«D dx ( 9 ) 

The plausibility of equation (9), at least for the limiting 
case of zero heat transfer, may be shown as follows:. Consider 
the flow in the boundary layer to be essentially two-dimensional, 
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inasmuch as the radius of curvature of the tube is large com- 
pared with the thickness of the boundary layer. Velocity and 
temperature gradients in the flow direction may be neglected. 
The change in the total energy of the fluid flowing through an 
elementary volume is the net heat conducted into the volume 
by the temperature gradient perpendicular to the flow direction 
plus the work done on the volume by the viscous shearing forces 
resulting from the velocity gradient perpendicular to the flow 
direction. In the form of an equation, the energy balance 
becomes 

c v ±61 , l(& + , 1 5 ] do) 

Sx By \ ay J dy 

If c , k, and \x are assumed to be independent of the tempera- 
ture, equation (10) may be rewritten as 

ST 




pVc = hc T + / (11) 



P 3x 



For the case t iC p ( /'^ ~ ^> a particular solution of equa- 
tion (11) is 

V 2 

T + =-= — = T = Constant ( 12 ) 

^jc p s 

The net heat conducted into the elementary volume is, for this 
solution, equal to the work done by the volume against viscous 
shearing forces, so that the total -energy distribution is 
uniform throughout the boundary layer. It is apparent that 

at the wall, where V = 0 and T = T > §*■ = 0: that is, the wall 

oy 

is at temperature T g and the heat transfer is zero. 

llC-p. 

The validity of the assumption £ = 1 depends upon whether 

k 

the boundary layer is laminar or turbulent. There will usu- 
ally be an outer turbulent boundary layer adjoining a laminar 
sublayer, part of the temperature rise occurring within each. 

For the laminar sublayer £ = 0.75. For the turbulent part of 

k 
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the boundary layer, however, — ~ m 1 because the effective values 

k 

of k and \i are greater for turbulent flow than for laminar ^ c 
flow (reference 4), The over-all effects of the variation of 

k 

and 



with the nature of the flow and of the variation of n, 



V 



k With temperature on the equilibrium temperature of the pipe wall 
as given by equation (11) are not very great. Frossel in refer- 
ence 5 observed that air accelerating from rest at atmospheric 
t€>mperature to supersonic velocities did not appreciably change 
the tube -wall temperature from atmospheric temperature , even 
though the temperature of the air dropped more than 100° F. 

In the development of the heat -transfer equation for high air- 
speeds, the kinetic -energy correction must be added to equa- 
tion (5) before this equation may be used in the high -velocity 
range; thus 



dQ = 



™ pV 



dT + d 



'6 



or 



dQ . SJ- P Vc p dT 



9 



(13) 



If the hsat -transfer coefficient is assumed to be the same at 
high speeds as at low speeds, the solution of equations (13) 
and (9) is analogous to equation (8): 



q = P V SE- c 
4 P 



t0 2 



1 - e 



.4(0.0247 )R 



-0.2, 



D 



(14) 



Equation (14) is recommended as a first approximation to the 
heat -transfer rate at high airspeeds for which experimental data 
are lacking. The assumption in the derivation that the heat- 
transfer coefficient remains the same at high speeds as at low 
speeds is implied in the rule given by Grocco in reference 4. It 
is an experimental fact (reference 5) that the skin-friction 
coefficient is independent of the Mach number; and, because the 
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mechanism causing heat transfer is essentially the same as the 
mechanism causing skin friction, it seems justifiable to assume 
in the absence of direct experimental data that the heat- 
transfer coefficient is also independent of the Mach number. 



Static-Pressure Losses in the Tubes 

Linear axial-vel oci ty distributions. - The primary cause 
of static -pressure drop in the flow through a hot radiator tube 
is the shearing force at the tube wall due to the skin fric- 
tion. According to reference 3, the static-pressure drop due 
to skin friction is 



_? m f) p Y 

Ap - 0.098pT J R 
F 



or, in the differential form with due regard for sign, 



~dp F = 0.098pV 2 ~ 



-0.2! 



dx 



(15) 



Because the product pV must be a constant within the 
tube, the density reduction within the tube results in an 
acceleration that causes the further pressure drop 



m dx 



D/2 



- 0 



J* (2jty) dy 



Ax 



When the variable y is changed to g^-> the expression 
rewritten in the differential form is 
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For the Blasi us —power-law velocity distribution the value of 
7 

the integral ia 50/49. A value of "1 will be assumed; that is, 
the average velocity is supposed equal to the root -mean -square 
velocity, in the re3t 6 the paper the bar has been omitted 
from V and no distinct: n Is indicated between velocity, 
average velocity, and root -mean -square velocity. 

The equation for the total loss of static pressure is given 
by the sum of equations (15) and (16) as follows: 



-dp = 0.098pV 6 



•0.2 



dx + 



(17) 



As pV is constant, 




(16) 



The static-pressure drop in the tube may be approximated if the 
axial-velocity distribution is assumed to be linear, an /assump- 
tion that will be justified in the following section. By this 
assumption 

By use of 



and 



P? T 2 P3 T 3 
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equation (19) may "be solved for po - p-z; in terms of the known 
quantities p«, Pg, V?' an< ^ L A>: 



P 2 " P3 = 



l (B -/b 2 - 4c) 



(20) 



where 



B 



?2 ^ ^2 -l)' 



C 




and 




F 



2 1 0.049 £ B 



The value of is obtained from equation (6). 

Nonlinear axial-velocity distributions, - When static- 
pressure drops exceed 30 percent of the absolute pressure at 
the tube entrance , an error of more than 10 percent will be 
made in determining the pressure drop from equation (20). In 
this case the simplified solution of the differential equation 
for the pressure drop based on the linear axial-velocity dis- 
tribution is no longer valid. A more precise solution of equa- 
tion (17) may be derived as follows: First } the axial velocity 
is expanded as a power series of the x/D ratio with coefficients 
to be determined. This expression is substituted in equation (17) 
and the equation is integrated for the pressure as a function 
of the x/D ratio. By the use of the series development for the 
stagnation temperature, the coefficients of the power series 
are finally determined . 

Expanding the axial velocity in terms of the x/D ratio gives 



V 9 + a| 



(!) ♦ iff ♦ 4f * <*) 4 • • • 



(21) 
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which' is substituted in equation (17) Rewritten as 



-dp = PV 



mV + 



dV 



where m = 0.098E" 0,2 . There is thus obtained for the axial- 
pressure distribution 



P = P 2 - pVm 



pV 



(22) 



If T is replaced by T, and 0 is defined as T w - T g , there 



is obtained from equation (6): 

-rti 



o = e „e 



From the definition of T 



i - a&U .JSJL • • 

A D / 2! 31 



(23) 



T w " T s ' T w - ^ " 2Jc 



and .'the gas law 



P p PB Q T 



8 may be written as 



9 = 1 



pV 



w R G pV 2Jc p 



(24) 
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Linear equations for the coefficients in the power series 
of equation (21) are obtained by substituting the expressions 
for V, p, and 9 given by equations =(21), (22), and (23), 
respectively, in equation (24) and setting the sum of the ooei 
fictenta of the various powers of the ~ ratio equal to zero: 

D 



pa - ml ™~ ■+ e 2 ) 




(3b . 



3mV 2 a a 2 (.7 + l) Sfcm 2 



2B, 



2 7 R r 




pd = 



ah (7 + 1) 9 2 m " 
— + 



7B 



G 



3! 



5raf V 2 c ab) y + X f9m/t . . 2s ?g£ 



7B, 



•ft- \ — f— + -T- I + ~s — (a* 3 + bc J + asr + -rt 



3 % 



V (25) 



where 



T. 



v 



■7B G 



The coefficients a, "b, c, d, and e are thus determined 
in terms of and the entrance conditions m, an( ^ ®2 % 

The -pressure drop for any value of x/D is found by substituting 
these values in equation (22) or in the following rearrangement 
of equation (22) : 



P 2 " P * P v 



(mV 2 + a)| 



ma 4- 2b 



mb + 3c 





(26) 
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The rapidity of convergence of the right-hand side of equa- 
tion (26) has "been investigated by a graphical integration of 
equation (17) for specific initial conditions. In figure 1 the 
results of this integration are compared with the results obtained 
by use of equation (26). Mathematical accuracy compatible with 
the other simplifying assumptions of this analysis is obtained 
for pressure ratios down to 0.55 by using five terms of the power 
series. 

Figure 1 is also useful for checking the accuracy of equation (20), 
which was derived by assuming a linear axial-velocity distribution. 
Equation (20) appears to be reasonably accurate throughout the present 
design range and down to a pressure ratio P3/P2 of °« 7 * 

Another approximate solution for equation (17), based on the 
use of the stagnation temperature in the gas law, is 



where 




The accuracy of this solution, as indicated in figure 1 is between 
the accuracies of the two solutions already discussed. 

Effect of Mach number on flow in radiator tubes. - Important 
information concerning the effect of Mach number on flow in radi- 
ator tubes can be obtained from equation (17). From equation (23) 



9 „ v -o.o98R-°'?j< ..• . '\ 

which, when differentiated, is 

= -0.098 1^! dx (28) 
6 D 
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Substituting directly in equation (17) results in 



-dp = (-pV)V ~ + pV dV (29) 

If dp is now eliminated from equation (29) by means of the 
two relationships 

P = P%T 

and 



T + J. = T 

2Jc p a 



the resulting equation in two variables is 



2JoJ +R g( 0 - T w) 



M. = S. \ p/ - : (30) 



Equation (30) is important because it is independent of geometry 
and Reynolds number. 

At present no closed solution of equation (50) is available. 
A graphical solution/ determined by the isoclinic method, is given 
in figure 2. The curve for any particular radiator is the curve 
that passes through the point ^ 2 2 ^ 2 ) corresponding to the 
given entrance conditions. The straight line in the figure sep- 
arates the supersonic range on the right from the subsonic range 
on the left. The equation of the straight line is found by 
setting d8/dV equal to zero in equation (30) 



(31) 
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This equation may be easily reduced, from the definition of 9, 

to - - 1; that is, the line corresponds to a Mach number of unity. 

8. 

If the entrance velocity V 9 is subsonic, V will increase 

d 0 

along the tube and may reach sonic velocity, where _ = 0. 

For the air to accelerate into the supersonic region,. 0 would 
have to increase; that is, the total energy would have to decrease, 
which is contrary to equation (23), the heat -transfer equation. A 
supersonic velocity may then never develop within a radiator tube 
with subsonic entrance velocity. The limiting Mach number of unity 
will be found at the tube exit. 

If the entrance velocity is supersonic, V will decrease 
along the tube, the change being now from right to left along the 
curves of figure 2. The flow may not, however, pass continuously 
into the subsonic region without violating equation (23). A 
discontinuous transition may occur from the supersonic to the sub- 
sonic region by means of a stationary compression shock if the 
upstream flow is supor3onic, depending on the tube length and the 
exit pressure. Such phenomena, however, are of no importance in 
current radiator technology and a discussion of them is beyond the 
scope of this paper. 

Pressur e lo33 at th e e.xitLJ^ " For radi ~ 

ators in current use the discharge of the air from the radiator tub 
is accompanied by a loss of total pressure of about O.Sq (ref- 
erence 2, p. 10). This loss somewhat exceeds that given by the 
well-known Borda-Carnot formula for expansion loss 



due partly to surface irregularities at the exit and partly to 
the difference in kinetic energies associated with the upstream 
and the downstream velocity distributions. 





is only 0.1 lq.^ 



The difference is probably 



O utlet flow. -' For convenience in computation, the air behind 
the radiator is assumed to be brought to stagnation conditions. 
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The stagnation pressure is then 



f 

p = p, 1 + 

3 4 



y t 1 3 V 1 = 



=0.2q, 

27 R G T 3 y > 



(33) 



and, from equation (6), the stagnation temperature is 



% - *s 2 * ( T W " T B 2 



-4(0.0247)B" 0,2 ^ 
1 - e r 



(34) 



The air is now assumed to expand isentropically to the free- 
stream pressure. The velocity of the air at this pressure is 

given by 



2Jc p T II - T 



s 4, 



(35) 



Introducing the relationship 



7-1 



T 



'V 



(36) 



allows equation (35) to be rewritten as 



P is / 



1 - 



7-1 



(37) 



The drag power per unit open area due to the momentum change 
of the cooling air is . 
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„ P2 V ^ (V o - V 5 )V o ' (38) 
P D = 550 

It is to "be noted that this power represents, in general, a prac- 
tical minimum as this value neglects any further effects of the 
low-energy cooling air .on the external drag. 



PERFORMANCE CHARTS 

For the radiator performance charts the pressure drops within 
the radiator tubes were, for convenience, calculated hy the simplest 
of the methods discussed, namely, the method that assumes a linear 
velocity distribution along the tube. As has already been indi- 
cated, the pressure losses are thereby slightly overestimated^ or 
those cases in which the pressure drop is a large fraction oi the 
absolute pressure. The calculation of the flow path of the cool- 
ing air to obtain the performance charts is outlined as lollows: 

Values of airplane speed, altitude, and tube-entrance velo- 
city are assumed. The temperature T 2 at the tube entrance 
follows from equation (l), but the pressure p g given by equa- 
tion (2) is arbitrarily reduced by 10 percent of the free-stream 
dynamic pressure to account for flow separation and skin friction. 
The heat-transfer rate follows from equation (14) by evaluating 
the viscosity in the Reynolds number at the mean air temperature 
approximated by use of equation (6). The pressure at the radiator 
exit follows from equation (20), where T 3 is approximated by 
equation (6). The stagnation pressure and the temperature in the 
duct behind the radiator follow from equations (33) and {oi), 
respectively. Finally, the drag power is obtained from equa- 
tion (38) with the exit velocity of equation (37). The atmos- 
pheric conditions were assumed to be those of Army air and the 
radiator -tube -wall temperature was assumed to be 240° F. A 
ratio of free -flow area to frontal area of 2/3 was assumed; 
for any other ratio, the results merely change in proportion. 

The calculations covered the following ranges of variables: 

... 0 to 50,000 

Altitude, ft 0 to 500 

Airspeed, mph . . . • 9 1? and 15 

Radiator-tube length, in • • • ■ > d Wg 

Radiator-tube diameter, in i/ • 
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The six performance charts of figure 3 are for the six 
different radiators. In these charts the total-pressure loss, 
defined as the difference between the stagnation pressures 
ahead of and "behind the radiator, lias been plotted against 
the heat -transfer rate per unit frontal area. The inlet pres- 
sure loss is not included. The ordinates have been divided by 
0 for NACA standard air (see table I) merely for convenience 
in separating the curves. Line3 of constant tube -entrance 
velocity hav? also been plotted in figure 3 to aid in the quick 
determination of the cooling-air quantity at the radiator 
face . 

Although the calculations were carried out for airplane 
speeds of 1.00 miles per hour to 500 miles per hour, the results 
plotted in figure 3 are only for an airplane speed of 300 miles 
per hour. The ratio of the heat -transfer rate at several air- 
plane speeds to the heat -transfer rate at an airplane speed of 
300 miles per hour for the same total -pressure loss has been 
tabulated for different altitudes in table II. These cor- 
rections, although appreciable in some cases, especially at the 
lower altitudes, are for the most part negligible. A similar 
correction for the cooling-air quantity is given In table III, 
where the ratios of the tube -entrance airspeeds at several 
airplane speeds to the tube -entrance airspeed at an airplane 
speed of 300 miles per hour for the same loss in total pres- 
sure are tabulated for different altitudes. 

It may be remarked that the corrections in table II involve 
two opposing effects of airplane speed: namely, increased tube- 
entrance stagnation temperature, which decreases the heat- 
transfer rate for a given total pressure loss, and increased tube- 
entrance density, which has the opposite effect. At low altitudes, 
the temperature effect predominates: but at high altitudes, where 
the critical design condition usually occurs, these effects 
almost compensate each other and render figure 3 particularly 
accurate . 

A comparison of figures 3(a) and 3(d) shows that about 50 per- 
cent more heat may be dissipated per unit flow area by -L-inch 

tubes 9 inches long than by L-inch tubes 9 inches long for the 

4 

same total -pressure loss. In general, for tube lengths of this 
order, the tubes of smaller diameter will permit smaller radi- 
ators and/or smaller total-pressure losses. 
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Figure 4 shows the drag characteristics of the radiators 
as functions of the flight speed, the heat -transfer rate, and 
the tube -entrance velocity. The heat -transfer rate and the tube- 
entrance velocity are obtained from figure 3. No Interference 
effects or weight drag are included. It may be seen that, at 
the high inlet velocities necessary at high altitudes with radi- 
ators of reasonable size, the cooling drag power becomes very high, 

In the past it was hoped that the Meredith phenomenon would 
bring about very low drags at high altitudes. Appreciable decreases 
in drag are due to this source; but the calculations, which auto- 
matically take into account the Meredith phenomenon, show that 
these decreases are not sufficient to keep the drag of small-size 
high-pressure -drop units from appreciably increasing at high 
altitudes. 

An example will illustrate the use of figures 3 and 4 and of 
table II. The airplane for which the radiator will be designed is 
assumed to have the performance shown in table I. The densities 
and Mach numbers U3ed'in calculating q c are for Army air. The 

radiator design is assumed to be for tubes of t-inch internal 
diameter and 9 -inch length. 

The usual condition determining the radiator frontal area 
will be either climb at sea level or climb at the maximum altitude, 
that is, the condition at which the heat-transfer rate per unit 
frontal area is lowest when the entire total pressure available 
for cooling is utilized. If the available pressure for cooling is 
assumed to be 0.9q , the values of the heat -transfer rate may be 

read from figure 3(a) for both high speed and climb. The results, 
corrected according to table II, are summarized in table IV. 

An examination of table IV discloses that the frontal area of 
the radiator is determined by the value of the heat -transfer rate 
in climb at 40,000 feet, if the necessary heat -transfer rate for 
satisfactory cooling is assumed constant. If 1000 horsepower is 
to be transferred to the cooling air, the necessary radiator 
frontal area is 1000/.164 or 6 square feet. 

Once the radiator frontal area is fixed, the operating line of 
the radiator is the ordinate in figure 3 through a heat -transfer 
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rate of 164 horsepower per square foot. The ratio; of the drag 
power to the heat -transfer rate may be determined from figure 4(a) 
from the values of the airplane speed, the tube -entrance 
velocity, and the altitude. The results arc tabulated in table V. 



CONCLUSIONS 

From the present study of heat transfer, pressure drop, 
and drag power of radiators in flight at high altitudes, it 
is concluded that: 

1. At high altitudes, the density reduction within the 
radiator tubes results in appreciable increases in both 
friction pressure drop and acceleration pressure drop. 

2. The usual heat-transfer equations may be retained at 
high Mach numbers, provided that stagnation temperature is used 
in place of actural temperature. 

3. Heat-transfer rate as a function of pressure loss is 
practically independent of airplane speed at high altitudes. 

4. Excessive drag is associated with the use of small 
radiators at high altitudes; slight increases in radiator size 
result in large decreases in drag. 

5. The Meredith phenomenon becomes insignificant when the 
radiator pressure drop approaches the total available pres- 
sure. 
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Table I. - assumed airplane performance 

Wing loading, 40 lb/sq ft; heat rejection, 1000 hp] 

I 



Altitude 
(ft) 


Relative 
density, o 
(a) 


^igh speed 


Climb 


Air- 
speed, V 
(mph) 


Free-stream 
impact 

pressure, q 
(lb/sq ft) 


0.9q c 
0 


Air- 
speed. V 
(mph) 


Free-stream 

impact 
pressure, q 
(lb/sq ft) 0 


0.9q 
0 


0 


1.000 


333 


275 


247 


192 


87 


78.3 


10,000 


.738 


367 


247 


301 


218 


84 


102.4 


20,000 


.533 


403 


220 


371 


246 


77 


130.2 


30,000 


.374 


450 


190 


456 


280 


70 


168.5 


40,000 


.245 


501 


170 


625 


354 


79 


290.0 



The values of o are based on NACA standard air as in fig. 3. 



TABLE II. - RATIO OF HEAT -TRANSFER RATES AT SEVERAL AIRPLANE SPEEDS 
TO HEAT -TRANSFER RATE AT AN AIRPLANE SPEED OF 300 MILES PER 



HOUR FOR IDENTICAL TOTAL-PRESSURE LOSSES, 7— f 

(Pc, 300 







^^Altitude 














v o ^) 


0 


10,000 


20,000 


30,000 


40,000 


50,000 


(mph) 














100 


1.08 


1.05 


1.03 


1.02 


1.00 


1.00 


200 


1.05 


1.03 


1.02 


1.01 


1.00 


1.00 


300 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


400 


• 93 


.95 


.97 


.98 


.99 


1.00 


500 


.83 


.88 


.92 


.95 


.98 


.99 



TABLE III. - RATIO OF TUBE-ENTRANCE AIRSPEEDS AT SEVERAL AIRPLANE 
SPEEDS TO TUBE-ENTRANCE AIRSPEED AT AN AIRPLANE SPEED 
OF 300 MILES PER HOUR FOR IDENTICAL TOTAL- 



PRESSURE LOSSES, 

(V 2) 3 00 





V( V 2> 30 o 


"\Altitude 
















0 


10,000 


20,000 


30,000 


40,000 


50,000 


(mph) 














ICO 


1.05 


1.05 


1.05 


1.06 


1.06 


1.06 


200 


1.03 


1.03 


1.03 


1.03 


1.04 


1.04 


300 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


400 


.96 


.96 


.96 


.95 


.95 


.95 


500 


.92 


.91 


.91 


.90 


.89 


.89 
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TABLE IV. - VARIATION WITH ALTITUDE OF MAXIMlj:,' POSSIBLE 



COOLING PER SQUARE FOOT 



Altitude 
(ft) 


Hiafa sneed 


Climb 


Airspeed 
(fror. table I) 
(mph) 


< P c>300 
(fron fig. 3(a)) 
(hp / sq ft) 




P c 

(hp / SQ ft) 


Airspeed 
(mph) 


(hp / sq ft) 


r c 


(hp / sq ft) 


(h c) 30 o 
(fron table II) 


(F c'300 
(from table II) 


0 


333 


305 


0.98 


299 


192 


190 


1.05 


199 


10,000 


367 


300 


.97 


291 


218 


200 


1.03 


206 


20,000 


403 


282 


.97 


274 


24G 


193 


1.01 


195 


30,000 


450 


235 


.97 


228 


280 


175 


1.00 


175 


40,000 


501 


195 


.98 


191 


354 


164 


1.00 


164 
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TABLE V. - PERFORMANCE OP RADIATOR OF 6-SQUARE-F00T FRONTAL 
AREA TRANSFERRING 1000 HORSEPOWER 





Vo 

(from fig. 3(a) ) 
(mph) 


High speed 


Climb 


Altitude 
(ft) 


v o 

(from table I) 
(mph) 


V p c 
(from fir*. 4(a)) 


P D 
(hp) 


(mph) 




P D 
(hp) 


0 

10,000 
20,000 
30,000 
40,000 


130 
135 
145 

172 

220 


333 
367 
403 
45C 
501 


0.102 
.072 
.064 
.104 
.176 


102 
72 
64 
104 
176 


192 
218 
246 

280 
354 


0.124 
.108 
.102 
.200 
.360 


124 
108 
102 

200 
360 



> 



r 
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Fig. 1 




Fig. 2 
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Fig. 3b 
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Fig. 3c 




Fig. 3f NACA ARR No. L4Illb 
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Fig. 4a 
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Fig. 4c 



Fig. 4d 
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Fig. 4e 




Fig. 4f 
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